Decoherence Dynamics in Low-Dimensional Cold Atom Interferometers 
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We report on a study of the dynamics of decoherence of a matter- wave interferometer, consisting 
of a pair of low-dimensional cold atom condensates at finite temperature. We identify two distinct 
regimes in the time dependence of the coherence factor of the interferometer: quantum and classical. 
Explicit analytical results are obtained in both regimes. In particular, in the two-dimensional (2D) 
case in the classical (long time) regime, we find that the dynamics of decoherence is universal, 
exhibiting a power-law decay with an exponent, proportional to the ratio of the temperature to 
the Kosterlitz-Thouless temperature of a single 2D condensate. In the one-dimensional (ID) case 
in the classical regime we find a universal nonanalytic time dependence of decoherence, which is a 
consequence of the nonhydrodynamic nature of damping in ID liquids. 
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Understanding quantum phases and phase transitions 
is a substantially more complex task than understanding 
their classical counterparts. In particular, in the quan- 
tum case equilibrium and nonequilibrium properties are 
inseparable and thus dynamical relaxation processes and 
dissipation may have a profound effect on the nature 
of equilibrium phases and phase transitions. Superfluid 
to insulator (SI) transition in disordered 2D supercon- 
ducting films or Josephson junction arrays (JJA) [l[ is 
probably the most prominent example. While it is of- 
ten thought of as one of the most basic and fundamental 
quantum phase transitions, many aspects of it are still 
very poorly understood. In this context it appears inter- 
esting to investigate nonequilibrium phenomena in low- 
dimensional superfluid systems in a simpler setting than 
is offered by such highly complex systems as amorphous 
superconducting films or JJA. Cold atom systems seem to 
be especially well-suited for such studies since they allow 
for an unprecedented control of the relevant experimen- 
tal parameters and for direct, real-time measurements of 
phase coherence [2, OJ 0] • 

In this Letter we present a study of the dynamics of 
decoherence in a system of two low-dimensional (ID or 
2D) cold atom condensates that are prepared from a sin- 
gle phase coherent condensate by splitting it, using an 
optical lattice or a radio-frequency-field-induced poten- 
tial on an atom chip. It has been recently demonstrated 
by a number of groups [1] that it is possible to split a 
single condensate in such a way that the phase coherence 
between the two halves is initially well preserved. This 
phase coherent state is, however, a highly nonequilibrium 
one for the split condensate system. The system will then 
relax to thermal equilibrium over time, in which the con- 
densates are completely incoherent and it is the dynamics 
of this decoherence process that we are interested in. 

The problem of decoherence in cold atom conden- 
sates has been studied theoretically before by several au- 
thors [5]. Previous studies mostly focused on the single- 
mode approximation, which assumes that only the lowest- 
energy phase mode needs to be taken into account. This 
approximation is reasonable in three-dimensional (3D) 



condensates, but breaks down in low-dimensional sys- 
tems, where one needs to take into account the whole 
continuum of hydrodynamic (low-energy) modes, respon- 
sible for the enhanced fluctuations [fj and the absence of 
conventional order. Quantum decoherence due to such 
continuum of modes has been explored in a recent work 
by Bistritzer and Altman Q. In this Letter we extend 
their analysis to investigate both quantum and thermal 
decoherence. We also discuss the role of the "squeezing 
factor" , determined by the finite splitting time. 

At a given temperature T we can divide the hydro- 
dynamic modes into two groups: classical, with energies 
e < T , and quantum, with e > T . The dynamics of each 
group of modes is governed by classical and quantum 
fluctuations correspondingly. As a consequence, we can 
also identify two distinct regimes in the time dependence 
of the decoherence process: quantum and classical. At 
short times, i.e. times, shorter than the inverse temper- 
ature t < 1/T, decoherence dynamics is dominated by 
the quantum modes and thermal fluctuations may be ne- 
glected. Since at short times the memory of the initial 
state of the split system is still preserved, the quantum 
decoherence process is strongly influenced by the nature 
of this initial state, which in turn is determined by the 
process, by which a single condensate is split into two. 
We adopt a simple model, in which during the splitting 
the system is assumed to be described by a Josephson 
Hamiltonian with a time-dependent tunneling term @,0]: 

H = f A r 2 - Ejco S (6) « ^Nl + ^-e 2 , (1) 

where N r — Ni — N2 is the relative number of atoms 
in the two condensates, which is conjugate to the rel- 
ative phase [8,N r ] = i. The charging energy E c = 
djJt/dN — n/N, where /i is the chemical potential, taken 
to be the same in both condensates, and Ej is the (time- 
dependent) tunneling energy. Eq.JT]) is the Hamiltonian 
of a harmonic oscillator with a characteristic frequency 
uj = \/E c E,j (we will be using Ti = ks = 1 units hence- 
forth). One can distinguish two stages of the splitting 
process While uj(t) > l/r s , where t s is the char- 
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acteristic time of the splitting, the splitting is approxi- 
mately adiabatic and the ground state of Eq. |T]) should be 
a good approximation for the actual state of the double- 
well system. Once w(t) < l/r s , the splitting process is 
no longer adiabatic. We can then approximate this sec- 
ond stage of the splitting process as instantaneous. This 
is equivalent to saying that we approximate the state of 
the fully split system, when Ej = 0, by the ground state 
of Eq.([T]) at the moment t* when the adiabaticity condi- 
tion fails, i.e. when uj(t*) = 1/t s . Then we obtain the 
following result for the relative number fluctuations of the 
initial state of the fully split system: AN r ~ yjN/fiT a , 
where N = gV = (N% + Nu)/^. This formula agrees 
well with the apparent value of AN r , which can be in- 
ferred from the measurements of Ref. The agree- 
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classical regime turns out to be universal, independent of 
the microscopic nature of damping and interactions. This 
can make our result especially useful for thermometry in 
2D condensates. 

We now provide the most important details of the 
derivation of the above results. We consider a system 
of two ID or 2D superfluids, which are prepared at time 
t = in a phase coherent state by, for example, split- 
ting a single superfluid in a double- well optical potential. 
This highly nonequilibrium state will then relax to ther- 
mal equilibrium. We will assume that this relaxation is 
entirely due to intrinsic processes in each of the two su- 
perfluids. We start from the following imaginary time 
action for the double- well system: 



ment becomes even better if we take into account the 
parabolic trap potential, which gives /i ~ TV 2 / 5 instead 
of /i ~ N for the box geometry. It is useful to repre- 
sent the above result in terms of a "squeezing factor" £, 
defined as l/£ = AN r /\/N — l/^//xr s . Our analysis is 
valid for £ 3> 1. When /it s < 1 one should take £ w 1. 
Our formula for AN r thus interpolates between the lim- 
its of fast splitting (r s < when AN r ~ V~N, and 
slow (adiabatic) splitting, in which case AiV r ~ 0. 

Choosing the initial state of the system to be a 
minimum uncertainty wavepacket of width l/AN r , we 
find that the coherence factor, defined as ^d{t) — 
(l/V)Re / dx (e ie ^ x '*^), where 6*(x, t) is the relative phase 
between the two condensates, V = L d is the d- 
dimensional volume of each condensate (d = 1,2) and 
the angular brackets denote both quantum and thermal 
averages, is given by: 
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Here a = 1,2 labels the two condensates (summation 
over repeated indices is implicit) and J is the resid- 
ual tunneling matrix element between the condensates. 
We will assume henceforth that J is negligible and set 
it to zero. Following Popov 13(, we rewrite the ac- 
tion using density-phase variables, which we define as 
*&er = > and expand about the uniform equal 

density saddle point q\ = Q2 = g = fi/g to obtain the 
following hydrodynamic imaginary time action: 
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where to ~ is a short-time cutoff. The Kosterlitz- 
Thouless temperature Tkt of a single 2D condensate and 
the Luttinger parameter K, in the ID case are defined 
explicitly below. In the fast-splitting regime (r s < 1/fi) 
in ID our result agrees with the result of Bistritzer and 
Altman 0. 

At long times t ^> l/T classical modes, with energies 
e < T, dominate the dynamics. In this case we find: 
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where to are (dimension-dependent) cutoff times (time 
dependences in Eq.(J5]) are valid when t > to), given ex- 
plicitly below. Two features are noteworthy here. In the 
ID case classical decoherence dynamics has a nonanalytic 
time dependence. As shown below, this is a consequence 
of the fact that in ID liquids damping at finite tempera- 
tures is always nonhydrodynamic, i.e. the damping rate 
is a nonanalytic function of momentum [lS^ |. In both ID 
and 2D cases the time dependence of decoherence in the 



where h a = g a — g. The saddle point expansion is le- 
gitimate at sufficiently low temperatures, such that the 
(exponentially small) contribution of phase slips (ID) or 
vortices (2D) to correlations can be neglected. The har- 
monic part of the action So describes undamped collec- 
tive density-phase modes of the double-well superfluid. 
The anharmonic part of the action Si is responsible for 
the dissipation and relaxation to equilibrium at finite 
temperatures, as will be shown below 13|, ll4|. As dis- 
cussed above, at short times t < l/T after the decoupling 
of the two condensates, the contribution of thermal fluc- 
tuations of the relative phase to decoherence dynamics is 
negligible and temperature can simply be set to zero: the 
dynamics of decoherence in this case is purely quantum. 
Furthermore, we may neglect the anharmonic terms in 
Eq.©, which in the quantum regime are unimportant. 
Introducing the relative and center-of-mass phase vari- 
ables 6 = 9i — 62, X — ($1 + $z)/2 and integrating over 
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density fluctuations, we obtain: 
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The relative and center-of-mass phase dynamics are thus 
completely decoupled at short times. As the action for 
the relative phase is purely harmonic, the problem of the 
relative phase dynamics at short times may be solved 
exactly 7]. Passing from the imaginary time action for 
the long- wavelength relative phase modes Eq. ^ back to 
the Hamiltonian, and rewriting it in Fourier space, we 
obtain: 
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where [6» k ,±l k /] = i<5 kk /. 

Taking the initial state of the split-condensate system 
|$o) to be a minimum uncertainty wavepacket, as dis- 
cussed above, it is straightforward to evaluate the time 
evolution of the coherence factor pj. One obtains: 
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and e k = y^gg/mk = c s k. Evaluating Eq.® at times 
t > 1/fX, we obtain Eq.([2]), with Tkt = Trg/2m and 
K, = (n/2)y^gjgm (these are weak-coupling expressions 
for Tkt and /C, which are expected to be accurate at low 
temperatures). Note that in Eq.([2]) we have explicitly 
separated the contribution of the k = mode (the com- 
mon factor), which will contribute separately from the 
k > continuum in a finite-size system. This is nothing 
but the "phase diffusion" which is the only contri- 
bution that exists at long times in bulk systems. 

We now extend the above theory to the classical, i.e. 
long time t > 1/T limit. The character of the long 
time dynamics is determined by the low-energy, classical 
modes, with energies e < T. We start from the hydrody- 
namic action for density and phase fluctuations Eq.(|5|). 
Unlike in the quantum case, the anharmonic terms in 
the action turn out to be crucial, as will become clear 
below. Despite the fact that the source of relaxation is 
purely intrinsic in our case, it turns out to be possible, 
and very useful, to cast the problem into the form [1 51 ] of 



an "observable macroscopic variable" interacting with a 
"thermal bath" of microscopic degress of freedom. This 
is possible thanks to the following observation. While in 
the imaginary time action Eq. ([5]) all degrees of freedom 
enter on a completely equal footing, the initial conditions 
for the relative and center-of-mass variables, introduced 
above, are drastically different. Indeed, we can formally 
model the process of splitting a single condensate into 
two by suddenly changing the value of the tunneling am- 
plitude in Eq. |(2J) from a large value at times t < to zero 
at t — 0. In hydrodynamic description this is equivalent 
to suddenly driving a large gap for the relative phase- 
density collective modes to zero. Center-of-mass modes 
however, feel this change only weakly, through anhar- 
monic terms in the hydrodynamic action. It thus seems 
reasonable to assume that the center-of-mass modes ap- 
proximately remain in thermal equilibrium throughout 
the splitting and subsequent relaxation process [16|. We 
can then think of these degrees of freedom as forming 
a thermal bath. The relative density and phase modes, 
in contrast, are strongly affected by the separation pro- 
cess and are far out of equilibrium at t = 0. We thus 
arrive at the picture of out-of-equilibrium relative phase 
collective modes, interacting with a thermal bath of the 
center-of-mass modes, the interaction being described by 
the anharmonic terms in Eq.©. To make this separa- 
tion of the degrees of freedom explicit we can perturba- 
tively integrate out the center-of-mass degrees of freedom 
(we expect this perturbation theory to work well at long 
times) in Eq.©, and obtain the following effective har- 
monic action for the relative phase variables only: 
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In 2D the lowest order (single phonon bubble) ap- 
proximation for the dissipative kernel II(k, ioS) is suffi- 



cient [lj, [17[. In ID the contribution of a single bubble 
diagram diverges "on-shell", i.e. at iu> = e k , which is a 
consequence of kinematics in ID, namely the fact that the 
laws of energy and momentum conservation are satisfied 
simultaneously [12]. It is then necessary to resum the 
most divergent (maximal number of bubbles) diagrams 
at each order [lj|. We obtain: 
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d= 2, 



(11) 



where 71 = \/2aTfC/irii (a sa 1.954 is a numerical con- 
stant) and 72 = 4:7rT 2 /3v3(j,Tkt are damping coeffi- 
cients, characterizing the strength of dissipation. Note 
that while 72 is dimensionless, 71 has dimensions of 1/ -y/e. 
The most interesting feature of Eq. (fTl"]) is the nonanalytic 
dependence of n(k, iio) on k in the ID case 12j, which 
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was first noted in a different context by Andreev. This 
means that damping in a ID liquid at finite temperatures 
is nonhydrodynamic, which is a consequence of the break- 
down of superfluid order in ID on length scales longer 
than the temperature-dependent correlation length [12j . 
It is this nonanalytic dependence of the damping kernel 
on momentum that leads to the nonanalytic time depen- 
dence of decoherence in Eq.((3]). 

Using Keldysh formalism, it is possible to show that 
the imaginary time action Eq. (110|) is exactly equivalent to 
the following real time quantum Langevin equation (lil . 
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where the quantum noise variable £(k, t) is defined 
by its autocorrelation function (C(k, fi)C( — k, £2)) = 
2II^(k, t\ — t%) and the functions Hrj are given by: 

IIi(k,w) = ilmR(k,iui — >w + iO+), 
IL R (k,oj) = i n 7 (k,w)coth(a;/2T). (13) 

Solving Eq . (fl2l with the initial conditions 9(k, 0) = 
0, 9(k, 0) = [2l|, by the Laplace transform, we obtain: 



(14) 
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Note that in the limit t — > 00 Eq. (fl4[) correctly repro- 
duces the equilibrium magnitude of thermodynamic fluc- 
tuations of 9(e). Now we can easily evaluate the time 
dependence of the coherence factor. At long times we 
obtain: 
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de v d {e)(9 2 (e,t)) d 



(16) 



where v d (e) is the phonon density of states in a d- 
dimensional superfluid. Evaluating the integral over e, 
we arrive at Eq.Q, with the cutoff times given by: 



f (3nfiK./T 2 , d=l, 
{ I/72T, d = 2, 



where = [8/(2a) 1 / 3 r(l/3)] 3 / 2 w 2.61. 
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